Polarized ensembles of random pure states 
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A new family of polarized ensembles of random pure states is presented. These ensembles are 
quite manageable and can be simply obtained from the unitarily invariant measure on the space of 
states by superposition of pure states. We will use the obtained results for two purposes: on the 
one hand we will be able to derive an efficient strategy for sampling states from isopurity manifolds. 
On the other, we will characterize the deviation of a pure quantum state from separability under 
the influence of noise. 
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I. INTRODUCTION 

In the last years many researchers have been investi- 
gating the typical properties of random pure states, i.e. 
unit vectors drawn at "random" from the Hilbert space 
associated to a quantum system. This subject has at- 
tracted the attention in several directions, and some im- 
portant results have been achieved mostly dealing with 
the characterization of entanglement [THTU] . 

The standard ensemble which has been intensively in- 
vestigated is that of random pure states with measure 
induced by the Haar measure on the unitary group. This 
ensemble, being the maximally symmetric one, imple- 
ments in a natural way the case of minimal knowledge 
on a quantum state [IT] . It is structureless, in the sense 
that the induced measure only depends on the dimen- 
sion of the total Hilbert space and it is not sensible to 
any tensor product structure [TJ [T2] . 

For this reason, the unitarily invariant ensemble is also 
known as the unbiased ensemble 10J . A natural question 
is whether this ensemble of pure random states can be 
used to construct different, more complicated, ones. 

Various approaches have been proposed independently 
by several groups that have introduced different phys- 
ically motivated measures on the space of pure states. 
Recently, Zyczkowski et al. j!3j have analyzed some kind 
of structured ensembles of random pure states on compos- 
ite systems that are invariant under local unitary trans- 
formations. From another perspective, De Pasquale et 
al. [2] have proposed a classical statistical mechanics 
approach in order to explore the isopurity manifolds of 
random states. In the same spirit, Mueller et al. |15) have 
recently investigated ensembles of random pure states 
with fixed expectation value of some observable, in the 
framework of the concentration of measure phenomenon. 
However, there are still many obstructions in carrying 
out these programs, and the links among them is not yet 
clear. 

This paper is intended as another step toward new sce- 
narios beyond the unbiased ensemble. This step is mo- 
tivated as an operational way to capture the isopurity 



manifolds, and turns out to be in particular cases sim- 
ilar to the structured ensembles proposed in [15] , Our 
idea is to move beyond the unbiased ensemble by using 
a natural operation at hand in the Hilbert space, namely 
superposition of vector states. In this work we will show 
that unitarily invariant measures interact nicely with the 
operation of superposition of states. On the other hand, 
as a remark, we want to stress that, in the large size 
limit, the robustness of Marcenko-Pastur's theorem [TB] 
prevents many of the potentially workable ensembles to 
deviate from the Marcenko-Pastur law. 

The paper is organized as follows. In Section III] we in- 
troduce the concept of polarized ensembles of pure states 
that will play a central role in our work. In particular, we 
will consider the deviation from the unbiased ensemble 
induced by the Haar measure on the unitary group acting 
on the whole Hilbert space of the system. In Section III 
we will characterize the polarized ensembles by using the 
local purity of a subsystem. This approach will be used 
both for the study of the ensembles and for the descrip- 
tion of an efficient procedure for sampling typical states 
from an isopurity manifold (and, therefore, with a fixed 
value of bipartite entanglement). In Section IV we will 



apply these results to the characterization of separability 
of quantum states under the influence of noise. Finally, 
in Section |V] we will draw some conclusions. 



II. POLARIZED ENSEMBLES OF PURE 

STATES 

Consider a bipartite quantum system described by a 
pure state in a finite-dimensional Hilbert space %. The 
bipartition in two subsystems (^4, B) will induce on % a 
tensor product structure such that 

H = Ha®H B - (2.1) 

We will consider, without loss of generality, the case with 

dimH A =N< dimU B =M, (2.2) 

whence dimT-L — NM. 



Let us focus on the general situation in which the 
(pure) state of the quantum system has the form 



M = a|Vi)+jShfe), 



(2.3) 



where \tpi) , \ipi) £ % are normalized random states, sam- 
pled according to arbitrary probability measures. Once 
the probability distributions of \tpi) and l^) are speci- 
fied, the random variable \ib) is characterized by a well 
defined distribution. 

Due to its privileged role in many context of random 
matrix theory [171 118] , in the following we will make large 
use of the unitarily invariant (unbiased) probability mea- 
sure on pure states, denoted by hmn, induced by the 
Haar measure on the unitary group U(H). Moreover, 
we will consider the product measure /xjv x Mm on pure 
states, which is left invariant under the action local uni- 
taries U^Ha) x U(Hb)- In the latter case, we recall that 
the space of pure states is foliated in orbits of local uni- 
taries labelled by the degeneracy of the Schmidt coeffi- 
cients. Each of these orbits is a natural domain for the 
locally invariant measure /ijv x Mm PS EQ] • 

In general, depending on the ensembles chosen for sam- 
pling |i/>i) and ^2), the state \ip) defined in Eq. (2.3) will 



exhibit very different properties. Let us briefly outline 
the relevant features in some cases of particular interest. 

If both \tpi) and \ib 2 ) are sampled according to the 
unbiased ensemble Hmn-, then \ib) is also a random state 
whose distribution is invariant under the action of the 
unitary group U(H) (independently on the values of a 
and (3). 

The opposite situation occurs when (^i) is a fixed pure 
state and \ip2) ~ (J-nm- In this case the weights a and 
/3 are relevant; if a 3> /3 the unbiased ensemble becomes 
"polarized" along the direction defined by |'0i)- This 
polarization phenomenon is of particular interest if one 
wants to study the deviation of the properties of an en- 
semble of quantum states from a fixed reference state. In 
particular, in this work we will specialize our analysis to 
the cases of |'0i) separable or maximally entangled (in 
the bipartite sense). 



III. TYPICAL LOCAL PURITY 



Local purity of one-parameter polarized 
ensembles 



In this section we will focus on the consequences of the 
polarization of the ensemble on the properties of bipartite 
entanglement between the subsystems A and B. As a 
measure of entanglement we will consider the local purity 
of one of the parties. Given a pure state \ip) £ H, it is 
defined by 



tab(IV')) = Tr^p 



PA 



Tr f 



(3.1) 



where pa is the reduced density matrix of party A. The 
upper bound ttab = 1 and the lower bound wab — 



1/ dim "Ha = 1/N correspond, respectively, to separable 
and maximally entangled states. 

Le t us consider a state \\b) obtained as a superposi- 
tion (2.3) where, in particular, a fixed pure state is su- 
perposed to an unbiased random one. We will get the 
following one-parameter ensemble 



= et AB + y/l - e 2 U A i 



(3.2) 



where the normalized state |0o) € "H is fixed, e £ [0, 1] is 
a tunable parameter, Uab is the identity operator, and 
Uab £lA {%) is a random unitary acting on H, sampled 
according to the Haar measure on the unitary group. The 
state 



Uab |0o) 



(3.3) 



is therefore a (unit) random vector distributed according 
to the unitarily invariant measure on pure states pnm- 

Notice that for e = one recovers the unbiased en- 
semble. On the other hand, values of e ^ play the 
role of an offset which parametrizes, as discussed in the 
previous section, the degree of polarization of the ensem- 
ble in the direction of |0o). We also observe that, given 
two independent and symmetrically distributed random 
states |0i) , \<j>2)i the expectation value of their overlap 
vanishes, i.e. 



E[( t h.\<h)]=0. 



(3.4) 



As a consequence, the normalization of the state \ip) 
in (3.2) is assured on the average, in the sense that 



E[<V#)]=E[|h/>|| 2 ] =1, 



(3.5) 



and deviations from the average are exponentially sup- 
pressed for large N, as shown in Appendix [A} 

We are interested in the conditional expectation value 
of the purity itab{\iP}) given a fixed state \<po) and a 
bias e. Due to concentration of measures, for large iV 
this quantity will be the typical purity of the polarized 
ensemble ( |3.2[ ). 

The density operator associated to the random pure 
state lib) reads 



{i>\ = e 2 |0o) (0o| + (1 - e 2 ) 
+ eVl-e 2 (l0o)(^| + k 



k\), (3.6) 



where \<f>) is given by (3.3). We will use the following 
notation: 



(T 



Tr B |0) (0| , 

Tr B \<b ) (0o I, 

Tr B (|0 O ) (01 + 10) (0o I) 



(3.7) 



By tracing over subsystem B and performing a straight- 
forward calculation, we obtain the purity (which is a ran- 



dom variable) 

ttab(IV')) 



T?aPa 



= e 4 Trcr 2 

+ e 2 (l 



o^l-e^Tra 2 
e 2 ) TrS o 2 

2 /l ,2\ 



2e 2 (1 - e 2 ) Tr(a of r) 
2eVl - ( 2 Tr(a S 0<f ,) 

2\3/2 



2e(l 



Trier Si 



0,9 J 



(3.8) 



B. Typical local purity 



We now evaluate the expectation value of the purity 
ttab- The computation can be easily done by making 
use of a Gaussian approximation. More precisely, we 
will consider random vectors \<f>) G 'H whose components 
in an arbitrary basis are independent complex random 
variables normally distributed, At (0, 1/NM) (the nor- 
malization of \<j)) is assured on the average). 

The Gaussian approximation is fully justified for our 
purposes. Indeed, in the large size limit, concentra- 
tion phenomena and the simultaneous convergence of the 
Gaussian measure to the unitarily invariant measure on 
the sphere [21 provide the typicality of our results (see 
Appendix [A| for further details). Thus, averages on the 
unitary group can be substituted with averages with re- 
spect to Gaussian measures. In this case, expectation 
values of any smooth quantity of interest /(|c/>)) can be 
easily estimated. 

We claim that the typical local purity of ( |3.2[ ) depends 
only on the local purity of the pure state \<pof, 



tto = n AB (\(j)o)), 



(3.9) 



and not on the particular vector |c/>o) with the given pu- 
rity. Indeed, a direct calculation with \<j>) a Gaussian vec- 



tor shows that the only non vanishing terms in Eq. (3.8) 
are 



E[Tra 2 ] = 

E[Tr(a a)] = 

E[Trcr 2 ] = 

E[TrS 2 ] = 



Tr a = TT , 


1 


N' 


M + N 



MN 



M 



(3.10) 
(3.11) 

(3.12) 

(3.13) 



By plugging Eqs. ( 3.10 1-( 3.13 1 into Eq. (3.8) we finally 



get the conditional expectation value of 7tab(|V')) 
tt ab = E[ir AB | |0o) i e] = e 4 7i"o + (l - e 4 ) MN ■ 

This is a central result of the paper. 

Notice that if we had performed the average using the 
unitarily invariant measure unm for \<j)), as in (3.3), the 



only difference with the above calculation would have 
been in the term 



EfTrcr 2 



M + N 

MN + 1 : 



(3.15) 



as computed by Lubkin I22| . The relative difference with 
the typical purity of the unbiased Gaussian ensemble 



Tunb = 7TAB e=0 



M + N 

MN ' 



(3.16) 



is thus of order 0(1/ (NM)) and negligible for large sys- 
tems. 



We point out an important consequence of Eq. (3.14) 



even if |<6 ) is substituted by a state \<j)' Q ) = Ua®Ub \<j>o), 
with Ua(b) <= U\Ha(B))i belonging to its local orbit 
(therefore evolving on an isopurity manifold with arbi- 
trary measure), the value of the typical purity given by 
Eq. ( 3.14[ ) is not affected. Incidentally, this can also be 
seen as a consequence of a fundamental property of con- 
ditional expectations: 



E[f(X,Y)]=E[E[f(X,Y)\Y}' 



(3.17) 



where / is a function of two random variables X and 
Y. We will make frequent use of this property in the 
following. 



C. Typical purity for a separable and a maximally 
entangled polarizing state \(j>o) 

In this section we will discuss two interesting examples 
of the behavior of ttab, given the value of ttq. In particu- 
lar, we will choose the extremal situations of a separable 
or maximally entangled polarizing state \4>o}- 

Let us start by considering the case of |0 O ) separable 
with respect to the bipartition (A,B), so that ttq = 1. 
Acco rding to the discussion at the end of the Section 
IIIB we can allow |</>o) to be a random pure separable 
state and consider the generalized ensemble 



m 



\h 



scp/ 



iU A ®U b + Vl - £ 2 U Ai 
\4>o) a ® \<Mb 



|0sc P ),(3.18) 
(3.19) 



where \4>o)a anc ^ \4>o)b are nxe d states in T-La and %b, 
respectively, Ua and Ub are random local unitaries and 
Uab is a random global unitary transformation. The 
typical purity (3.14) of the one-parameter ensemble with 



a separable polarization (3.19) reads 



EIttab 



M + N 



+ f 



MN -M-N 



-sop/ , cj MN , . MN 

(3.20) 
The other extreme case is given by \4>o) = |</Wt), a 
maximally entangled pure state, so that ttq = 1/N. Then 
we are dealing with a polarized ensemble of the form 



eU A <Z)U B + \/\ - e 2 U Ai 



|</>ent) 



(3.21) 



where the reference state is such that 



Tr B |</)ent) (</>ent| = "at^A- 



(3.22) 



Such a polarization decreases the typical purity of the 
unbiased ensemble to the value 



E[tTAB I |^ent) > e] = 



M + AT e 4 



MTV M 



(3.23) 



These results will be compared in the following section 
to a numerical Montecarlo approach. 



D. Generation of random pure states with fixed 
local purity 

The results obtained in the previous sections suggest 
a very inexpensive strategy for generating random pure 
states with fixed value of the purity ttab- Indeed, the 
numerical sampling of the uniform measure on the man- 
ifold of pure states \ip) parametrized by a fixed value of 
itab will proceed trough the following steps: 



1. Choose e e [0,1] such that 



KAB 



4 

e 7T 



1 



7t"unb i 



(3.24) 



where ttq = 1 or 7r = 1/JV if the desired value 
of itab is, respectively, larger or smaller than the 
unbiased typical value 7r un b in (3.16). 



2. Generate a pure state \ip) by superposition 



m = 



yr^ 



(3.25) 



where \cj>) is sampled according to the unbiased 
measure /Ujvm and \<fio) is a separable or maxi- 
mally entangled pure state sampled randomly ac- 
cording to the invariant measure under local uni- 
taries Hn><JIm j depending on the value of ir chosen 
in Eq. ( |3.24[ ). 



In Figure II] the analytical formulas (3.20) and (3.23) 
are compared to the Montecarlo results for the values 
of itab obtained by sampling pure states through the 
procedure outlined above. The comparison shows clearly 
the efficiency of the sampling procedure in providing the 
correct behavior of the typical purity vs the bias e with 
quite small fluctuations already for dimensions N — M — 
20. 



IV. ROBUSTNESS OF SEPARABILITY UNDER 
RANDOM PERTURBATIONS 

As an application, in this section we will use the re- 
sults obtained from the study of polarized ensembles to 
analyze the stability of "separability" of quantum states 
with respect to random additive perturbations. More 
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FIG. 1: (Color online) Typical purity of polarized ensembles 
vs e 4 , depending on the chosen reference state \ipo). We com- 
pare the analytical prediction (continuous lines) with the val- 
ues of ix ab obtained form the sampling procedure described 
in the text. We have considered balanced bipartitions with 
N — M — 20. Right side: the continuous line represents the 
analytical prediction for the purity of an ensemble polarized 
by a separable state [Eq. ( |3.20 l]. Depending on the value of 
the parameter e the purity ranges from the unbiased value 
7T un b = (M + N)/MN to the maximum tto = 1. Left side: 
the continuous line represents the analytical prediction for 
the purity of an ensemble polarized by a maximally entan- 
gled state [Eq. ( J3.23 1] . As e increases the typical value of the 
purity decreases from 7r un b to the minimum tvq = 1/JV. The 
error bars represent the standard deviations of the numerical 
simulations from the estimated average. Such fluctuations are 
exponentially suppressed as the dimensions N, M increase, as 
shown in appendix [A] 



precisely, if the state of the system \<j> ) e Ha <S> "Hb 
is separable, \4>q) = \(,a) A ® |xo)s, how much noise is 
necessary to make the reduced state pa = Trs[|0 o ) (<£o|] 
distinguishable from a pure state? We notice that the 
problem of the characterization of separability has at- 
tracted a lot of interest in the context of the analysis of 
ground states of quantum spin systems [2"3Tf2"S"] . 

Let us consider the state of the bipartite system in the 
form 



v/T 



V 2 \$o) a ®\xo)b+V I 



< 77 < 1, (4.1) 



where \<f>) ~ (J-nm is an unbiased random perturbation, 
and 77 measures the strength of the noise. (Notice that 
in the language of polarized ensembles of the previous 
sections we have e = \J\ — r/ 2 ). 

This setup models a general situation: some non- 
controllable noise prevents one to prepare with infinite 
precision a given state. The noise has the consequence 
that the reduced density matrix pa is not a projection but 
is mixed. It is reasonable that, for weak noise (77 ^C 1), 
we still obtain a reduced state pa that is close to a pro- 
jection. Then the question is: is there a threshold on the 
value of noise above which separability is appreciably de- 
stroyed? 



In order to answer this question, it is convenient to 
introduce the notion of effective dimension (also known 
as "participation ratio" [55]) of a state p defined as 



d eS (p) 



1 



TV 



e[i,N] 



(4.2) 



The effective dimension of a mixed state quantifies how 
many pure states appreciably contribute to the mixture. 
Moreover, differently from the rank of p, d cS captures 
the probabilistic weight of different states and is more 
manageable for explicit calculations. 

For a separable pure state, \£o) A ® IXo)b> the reduced 
density matrix, pa = |£o) (£o|) nas effective dimension 
d ■ (/9a) = rank(p J 4) = 1. A global perturbation acting 
on |£o)a® \xo)a can ^ e appreciated locally if the reduced 
state pa becomes a mixture. In order to obtain a mixture 
one needs at least two pure components. Therefore, we 
can say that a quantum state is distinguishable from a 
one-dimensional projection if its effective dimension d 
is equal or larger than 2. As a consequence, we are led 
to the following criterion on the separability of the state 
averaged over the noise realization: 



V. CONCLUSIONS 

In this paper we have shown that, using the superpo- 
sition principle, we can take advantage of the knowledge 
of the unbiased ensemble of random pure states in or- 
der to explore new interesting ensembles. In particular, 
we have found that adding a bias in a suitable direction 
is enough to polarize the Haar invariant ensemble. We 
stress that our approach has been oriented to the study 
of typical bipartite entanglement between subsystems, as 
measured by the local purity. 

This strategy yields an efficient and simple sampling of 
random pure states with fixed value of purity, and paves 
the way to further explorations and a deeper characteri- 
zation of the geometry of isopurity manifolds. 

Finally, we have applied our results to the analysis of 
separability of quantum states under the influence of ran- 
dom perturbation modeled through a coherent superpo- 
sition. For large systems, we have obtained a critical 
value of the noise strength, independent of the system 
size, beyond which the state is no longer separable, and 
the reduced state gets appreciably mixed. 



d ctt (E[p A })<2. 



(4.3) 



We get that 
d eS (E[p A }) = 



E[Trp A _ 



E[ttab| |0 SCP ) , y/l - V 2 ] ' 

(4.4) 
and then, from (3.20), it is straightforward to prove that 



condition (4.3 1 is satisfied if 



V 2 < ?7*(7Tunb) = 1 



1 - 27T unb 

2 - 27T unb ' 



(4.5) 



where 7T un b is the typical purity of the unbiased ensem- 



ble (3.16). If r\ > 77* the local state pa will be mixed with 



high probability. 

In the limit of large system sizes, N, M 
threshold critical value becomes 



r)l 



1 

72 



o 



-co, the 



(4.6) 



since 1/iV < 7r un b < 2/N. Therefore, as long as the state 
j^) of the large quantum system has the form (4.1 ) with 



T}< 



0.54, 



(4.7) 
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Appendix A: Gaussian approximation and typicality 
of the average purity 

The typicality of the average purity in the polarized 
ensembles used in Section |IIIB| relies on the following 
concentration phenomenon for Gaussian variables |I8j : 

Lemma A. 1 . Let X = ( X\ , X2 , ■ ■ ■ , Xk ) be a vector with 
independent identically distributed Gaussian components, 
with distribution Xi ~ A/"(0, a 2 ). Then, for any smooth 
function f : K fc — > R, with r\ — sup |V/| < oo, the follow- 
ing concentration inequality holds 



Pr{|/(X)-E[/]|>a}<2exp 



a 



4?7 2 cr" 



(Al) 



Le t us now consider the polarized ensemble defined in 
Eq. ((O): 



one has d oS (pA) < 2, and separability will be (approxi- 
mately) preserved. Notice that in this case, by applying 
perturbation theory, one gets that the spectrum of pa 
is made of a large eigenvalue of order 0(1 — rj 2 ) and a 
sea of eigenvalues of order 0(1/N) that have a negligible 
influence on the reduced density matrix pa- 



LAB 



+ Vl- e 2 U Ai 



< e < I . (A2) 



In the Gaussian approximation the 2MN real coordi- 
nates of the random vector state are Gaussian i.i.d. ran- 
dom variables with distribution Af (0, 1/2M N) (so that 
normalization is assured on average). The local purity 



function ttab, denned in (3.1), has Lipschitz constant 



bounded by rj < 4, and then by applying the Lemma one 
easily obtains 



Pr\ 



| TT AB {\iP))-E[Tr AB } >a|<2exp 



NMa 2 



32(1 



(A3) 

Incidentally, we mention that a similar Gaussian tail 
can be derived for uniformly distributed unit vectors by 
Levy's lemma [57]. The proof relies on a judicious use 
of (5-nets. In a finite dimensional setting, any subset of 
the sphere of states is totally bounded, in the sense that 



it admits a finite 6- net, for all 6 > 0. What we are 
interested in is a bound on the cardinality of a S-net N 
on manifolds of equal Schmidt rank k. For such manifolds 
a bound is given by [55] 



W\< 



2k(N+M) 



(A4) 



This estimate is good enough to control the probability 
of deviations of the purity from its average, so that a 
bound of the form (A3) is obtained. 
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